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1 Introduction 



An oft-repeated statement is that the recent evidence for neutrino mass and oscillations jT] 
constitutes proof of physics beyond the Standard Model. The important sub-text to this, 
however, is that the nature of this physics remains, at present, largely unknown, with very 
many possibilities and open questions. The quest for this new physics is closely linked to 
several unknowns in the neutrino sector : 

(a) The absolute scale of neutrino masses, 

(b) The hierarchy of neutrino masses, 

(c) The magnitude of the lepton mixing matrix element Ues (see below), 

(d) The value of the leptonic CP violating phase 5cp, and 

(e) The Dirac or Majorana nature of neutrinos. 

The hierarchy of neutrino masses rrii, ((b) above) with i = 1,2,3, refers to the ordering 
of the mass eigenstates. Given the fact that A21 = — ml > from solar neutrino data, 
knowing the hierarchy translates to determining the (as yet unknown^) sign of A31 = m| — m^, 
henceforth referred to as sign(A3i). Specifically, if this sign is positive, {i.e. m\ > > ml) 
the hierarchy is termed normal (NH), and if it is negative {i. e. m\ > m\ > m|), the hierarchy is 
said to be inverted (IH). The neutrino mixing matrix U, with elements Uai relating the flavour 
or weak interaction eigenstates (labeled by a = e, /i, r) to the mass eigenstates (labeled by 
i = 1,2,3) leads to the same flavour composition for a given mass state whether the hierarchy 
is normal or inverted. Attempts to construct a unified theory beyond the Standard Model 
however, depend crucially on sign(A3i); in fact, one way to classify families of models is via 
the neutrino hierarchy they assume as input. 

A large class of Grand Unified Theories (GUTs) use the Type I seesaw mechanism 
requiring the existence of heavy right-handed neutrinos at the GUT scale to generate small 
neutrino masses. Since GUTs typically unify quarks and leptons, and the quark hierarchy is 
normal, such models favour a normal neutrino hierarchy. Typically, if one uses an inverted 
hierarchy to construct a Type I seesaw model, one finds that the light neutrino mass 
spectrum may become very sensitive to small changes in the heavy neutrino masses and to 
radiative corrections. Inverted hierarchies imply the near degeneracy of the states m2 and 
mi, and since a corresponding degeneracy is absent in the quark sector, they may require 
the presence of an additional global symmetry in the lepton sector. Additionally, important 
features of leptogenesis, which may help partially understand the present matter-antimatter 
asymmetry of the universe, are significantly different in the two cases. In general, GUTs based 
on Type I seesaw models strongly depend on a normal mass hierarchy to obtain their many 
desirable features jHj. An inverted hierarchy, on the other hand, would favour the possibility of 
a unified theory based on a Type II seesaw mechanism, employing additional Higgs triplets, 

^Our knowledge of A31 — m| — derives from atmospheric and long baseline accelerator (K2K) data, 
which are sensitive only to its magnitude. 
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or on models which require a (global) lepton flavour symmetry. The hierarchy is thus a crucial 
marker in our quest for a unified theory, and its determination would be helpful in eliminating 
or at the very least strongly disfavouring large classes of such theories and in considerably 
narrowing the focus of this search. 

Generally speaking, determination of the mass hierarchy requires the observation of res- 
onant matter effects {i.e. long baselines) and a not too small sin^2^i3( ^ 0.05)^. Among 
the next generation long baseline accelerator experiments, combined results from T2K [5, ^ 
and the NuMI Off-axis experiment NOuA [7j may be able to infer the neutrino mass hierarchy 
|Hl m Uni m] • it is also possible to determine the hierarchy in experiments using beta beams 
^2] and neutrino factories fH]- However these experiments use the z/^ ^ z/g channel, the 
sensitivity of which is compromised by the ambiguities resulting from inherent degeneracies. 
Specifically, these originate in inter-relations between the sign(A3i), the phase 6cp and the 
mixing angle di^ [2 . To overcome this, the synergistic use of two experiments, sometimes 
with more than one measurement each will be necessary ITHl [T71 ITHl IT^ 1201 ■ Alternatively, 
two detectors at different distances (211 122] have been suggested. For baselines contemplated 
for the next generation of accelerator and superbeam experiments {i.e. ^ 800 km), much of the 
decrease in sensitivity arises from the Sep - sign(A3i) degeneracy. Its effect tends to decrease 
for much longer baselines, vanishing, in fact, for the magic baseline where the 6cp terms go 
to zero 123 123 • The use of earth matter effects in atmospheric muon neutrinos for the 
determination of the hierarchy has been studied in the context of magnetized iron calorimeter 
detectors in fl^ 123 EH] and water Cerenkov detectors in [20! ■ In particular, it was shown 
in [2H1 that the effect of degeneracies are not significant in the muon survival probability for 
the large baselines involved. Matter effects in supernova neutrinos can also, in principle, de- 
termine the neutrino mass hierarchy However large uncertainties in supernova neutrino 
fluxes reduce the sensitivity. Finally, if neutrinos are Majorana particles then determination 
of mass hierarchy may be possible from next generation neutrino-less double beta experiments 
provided hadronic uncertainties in the nuclear matrix elements are reduced from their 
present levels. Non-oscillation probes of the neutrino mass hierarchy and the limit of small 
sin^ 26*13 have been discussed recently in j^ . 

Megaton water Cerenkov detectors constitute an important future class of detectors with 
significant capabilities. The projects under active consideration are UNO jHHj in the US, a 
megaton detector in the Frejus laboratory fS^ in Europe, and the HyperKamiokande (HK) 
project inilSSl in Japan. Such massive detectors would enable impressive measurements of at- 
mospheric neutrino parameters. Specifically, they would determine | A31I and sin^ 623 to the few 
percent level, and obtain a much improved bound on ^^13. Since the water Cerenkov technique 
is not sensitive to the charge of the produced lepton on an event by event basis, the combined 
signal for neutrinos and antineutrinos must be searched for using statistical discriminators 
in order to determine the mass hierarchy. At present, the method which is contemplated for 
discriminating between neutrino and antineutrino interactions involves using the differences 
in the total cross section cr and its rapidity^ dependence {da/dy). Neutrino- nucleon interac- 
tions in the few GeV to 10 GeV range have a higher average rapidity than the corresponding 

^While these conditions are sufficient, they are not necessary. If they are not satisfied, however, quahtatively 
different experimental approaches entailing a significantly higher degree of precision will be required @]. 
■^Rapidity is defined as y = (E^ — Eicpton)/Ei/. 
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antineutrino-nucleon interactions, and thus produce more multi-hadron events. Thus an en- 
hanced signal for multi-ring electron-like events is expected in water Cerenkov detectors if the 
hierarchy is normal than if it is inverted [211 ES] • Given the importance of the hierarchy to 
theoretical efforts towards unification, it is recommended that we explore multiple approaches 
to this problem. In this paper, we exploit the sensitivity of the muon survival rate to resonant 
matter effects (as opposed to merely matter enhanced effects) for detecting the hierarchy in 
megaton Cerenkov detectors. Our method consists of a careful selection of energy and zenith 
angle bins for which the resonant matter effect in muon neutrino survival probability is ob- 
servationally significant. We perform realistic event-rate calculations for atmospheric muon 
rates in such detectors, which incorporate their efficiencies and lie within their resolutions. We 
show that for sin^ 26'i3 = 0.1 and ~ a 3 year run, signals with statistical significance in excess 
of 4(j are possible for pathlength and energy ranges where the muon survival probability is 
matter sensitive. As mentioned above, the results are largely free of parameter degeneracies 
and associated ambiguities. 



2 Earth - matter effects in P^^ 

As is well known the matter effects in neutrino oscillations arise due to the difference in the 
interactions of v^/vt and z/g as they traverse matter [32]. The mixing angle ^13 parameterizes 
the admixture of z/g with v^/vt in the oscillations driven by the larger mass-square difference 
A31. Hence all matter effects in oscillations, such as those which manifest themselves in long 
baseline and atmospheric neutrino experiments, are proportional to sine of this angle. The 
CHOOZ experiment constrains sin^ 26'i3 to be < 0.2 [^. For the energy which satisfies the 
resonance condition [3H] 

A31 cos 2^13 = 0.76 X 10-V(gm/cc) E(GeV) (1) 

the matter dependent ^13 gets amplified to 7r/4 and matter effects lead to large changes in the 
survival/oscillation probabilities, provided the pathlength dependent oscillating term is also 
close to 1. 

Resonant matter effects involving transitions have been extensively studied in 

the literature ^j. Most studies of next generation accelerator neutrino experiments, seeking 
to probe the mass hierarchy seek to exploit the matter effects in P^e- These experiments 
have baselines ranging from 300 km to 3000 km and energies in the ~ GeV range. For such 
'short' baselines and for energies in the few GeV range, while the percentage change in P^e 
can be large, its absolute value remains modestly small, since resonant matter effects do not 
develop. However, for the baseline range 6000 - 10500 km, since sin^ 29"^ then gradually builds 
to ~ 1, P^^ and P^r exhibit large matter sensitive variations, expecially if resonance occurs 
in the neighbourhood of a vacuum peak. In such a situation, the oscillating factors in the 
survival/oscillation probabilities are large. This in turn synergistically converts a large change 
in sin^2^^i3 due to matter effects into a large change in survival/oscillation probabilities. 

It is difficult to observe the resonant amplification of the matter effects in z/^ — > oscilla- 
tions, because the pathlength for which this occurs is inversely proportional to tan 2^13. For 
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almost all allowed values of 6^13, except those near the current upper bound, this pathlength 
is close to or larger than the diameter of earth jlUj : 

where p is the average density of the earth along the path and p is a positive integer > 0. 
However, the pathlength for which resonant matter effects in P^^ are maximum is given by 

HO] 

[pL]™;^ ~ pir cos 2^13 X 10^ km gm/cc, (3) 

where p is a positive integer > 0. For all allowed values of 6'i3, cos 2^13 is close to 1. Thus the 
resonant amplification in P^^ occurs within a narrow range of pathlengths for all the allowed 
values of ^13, thus making P^^ a suitable observable to study it. We obtain this pathlength 
to be about 7000 km by substituting the mantle density of earth p = 4.5 gm/cc and p = 1 
in Eq.Q. From a numerical study we have identified the energy range 5-10 GeV and the 
pathlength range 6000 - 9000 km to be suitable for studying resonant amplification of matter 
effects in P^^. In the case of P^r, the pathlength for which matter effects cause the largest 
change is given by jlHl ^ 

[pL]p" ~ I^^^-LUtt cos 2^13 X 10^ km gm/cc. (4) 

Here again the resonance amplification occurs for essentially the same pathlength for all al- 
lowed values of 6^13. Setting p = 1, in Eq.(jH), we get this pathlength to be about 9700 km^. 
A numerical study shows that the energy range 4-8 GeV and pathlength range 8000 - 10500 
GeV is suitable to study matter effects in P^^ induced by large matter-induced changes in 
p 

For the long pathlengths under consideration here, we need to explicitly take into account 
the varying density profile of the earth. We use the density profile given in Preliminary 
Reference Earth Model (PREM) |3T] . In Figure Q we plot P^^ as a function of energy for 
four different pathlengths in the range 6000 to 9000 km. In all of these four cases, the most 
significant matter effects occur in the energy range 5-10 GeV. In Figure El we plot P^^ as a 
function of energy for the two pathlengths 9700 km and 10500 km. In both these cases the 
most significant matter effects occur in the energy range 4-6 GeV. The curves Figure ^ and 
Figure El have been obtained by numerically solving the full three flavour neutrino propagation 
equation through earth matter. In obtaining these curves we have used the following values for 
neutrino parameters: IA31I = 0.002 eV^, A21 = 8.3 x 10"^ eV^, sin^^ia = 0.28, sin^^as = 0.5 
and sin^ 2^13 = 0.1. For each of the six pathlengths, the muon neutrino survival probability 
P^^ is calculated in vacuum and in matter for both the signs of A31. 

The following comments are in order : 

• For negative A31 (or inverted hierarchy), there is no discernible difference between vac- 
uum and matter survival probabilities. 

■*Note that in 001 cos 2^13 is misprinted as cos^ ^13. hep-ph/0408361, v3 has the correct expression. 
For p = 0, or the tt/2 peak, the matter-induced change in P^,- is seen to be small |2H]. 
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Fi gure 1: The muon survival probability P^^ in matter plotted as a function of E (GeV) for four pathlengths, 
L — 6000, 7000, 8000 and 9000 km and for the two signs of A31. The values of parameters used are sin^ Qy2 = 
0.28, sin^ 6*23 = 0.5, sin^ 26*13 = 0.1, IA31I = 2 x 10"^ and A21 = 8.3 x lO^'^ eV^. The survival probability 
in vacuum is also shown for comparison. 



• For positive A31 (or normal hierarchy), the value of P^^ in matter shown in Figure Q 
suffers a drop with respect to its vacuum value in the energy range 5-10 GeV. Whereas, 
the P^^ shown in Figure |21 undergoes an increase in the energy range 4-7 GeV. 

• For the pathlength range 6000 - 9000 km, the change in P^^, due to matter effects, is 
dominated by change in P^e- In the vicinity of the vacuum peak, matter P^^ is smaller 
by about 40% compared to the vacuum value. Three fourths of this change occurs due 
to change in P^e and the rest is due to change in P^^ • This is illustrated in Figure ^ 

• For pathlengths ^ 9000 km the matter effect in P^^ also becomes significant. For such 
pathlengths, there is a drop in P^,- which is as high as 70%, in the energy range 4 - 
6 GeV. This drop in P^^ overcomes the rise in P^e- Thus the net change in P^^ is an 
increase of the matter value over its vacuum value piHl28j. This is illustrated in FigureEl 
for two pathlengths 9700 km and 10500 km. 

• Beyond 10500 km the neutrinos start traversing the core and the mantle-core interference 
effects set in ||42j. It was shown in |l2Sj that at such pathlengths and energies relevant 
for atmospheric neutrinos the difference between the vacuum and matter event rates are 
not significant^. Hence we do not include such pathlengths in our analysis. 

^In the energy range ~ 3-6 GeV there is some difference between matter and vacuum rates due to mantle- 



6 



5 10 15 5 10 15 




5 10 15 5 10 15 



E (GeV) 

Fi gure 2: The muon survival probability P^^ in matter plotted as a function of E (GeV) for two pathlengths, 
L — 9700 and 10500 km and for the two signs of A31. The values of parameters used are sin^ dyi ~ 0.28, 
sin^fiiza = 0.5, sin^ 26*13 = 0.1, IA31I = 2 x IQ-^ eV^ and A21 = 8.3 x IQ-^^ eV^. The survival probability in 
vacuum is also shown for comparison. 



• In the case of muon antineutrinos, P^^ is essentially unchanged for the case of positive 
A31. It will experience resonance enhanced suppression for the case of negative A31, 
which again can be as large as 40% for pathlengths in the range 6000 - 9000 km. 

It is clear from Figure ^ and Figure |21 that matter effects in P^^ can be large, and a high 
statistics measurement of the muon survival rate in the energy and pathlength ranges indicated 
can be used to detect their presence as well as pin down the mass hierarchy. Atmospheric 
neutrinos passing through earth's mantle have the relevant pathlengths and possess energies 
in the desired range and thus are well-suited for this purpose. 

In general, these effects are sensitive to sin^ ^13, as mentioned earlier. This sensitivity, 
along with other aspects, was recently studied for a charge discriminating detector in j^El 
1^ . In [2H1, it was emphasized that a conclusive and statistically significant determination of 
the hierarchy and associated matter effects in such a detector requires a careful selection of 
pathlength and energy range. It was shown that a 4(T signal for matter effects is possible with 
an exposure of 1000 Ktyr in a typical iron calorimeter detector |S1 when muon events in 
the energy range 5-10 GeV and the pathlength range 6000 - 9700 km are considered. Wider 
ranges of energies and pathlengths result in an averaging out of the signal for matter effects 
pS] . In what follows, we study the sensitivity of the v^^ + z/^ survival rate to the hierarchy 

core interference effects |43| . 



7 



and to matter effects for a megaton size water Cerenkov detector. 

3 Numerical results for water Cerenkov detectors 

In order to determine the type of neutrino mass hierarchy, we exploit the muon neutrino and 
antineutrino disappearance channels. For a normal hierarchy (NH), the resonant amplification 
of matter effects occurs only for neutrinos, while for an inverted hierarchy (IH) it occurs only 
for anti-neutrinos. From Figure Hand Figure El we see that, for a NH, the /i~ rate undergoes a 
change due to matter effects but the /i"*" rate is essentially the same as the vacuum oscillation 
rate. In the case of an IH, the situation is reversed. Here we have to make a distinction 
between two different kinematic choices: 

1. E = 5 - 10 GeV and L = 6000 - 9000 km : In this energy and pathlength ranges the de- 
crease in P^^ is induced mainly by the increase in P^e- For a NH, matter effects suppress 
the /i~ event rates by a large fraction (about 25%) compared to their vacuum oscillation 
rates whereas the /i"*" event rates will be essentially the same as their vacuum oscillation 
values. For an IH, the /i~ rates are equal to their vacuum rates and the /i"*" rates undergo 
significant suppression. These results follow from the plots of P^^ shown in Figure H 

2. E = 4 - 6 GeV and L = 9000 - 10500 km : In this energy and pathlength ranges P^^ in- 
creases due to a sharp fall in P^,-. For a NH, matter effects increase the fi~ rates 
compared to their vacuum oscillation rates whereas, once again, the fi^ rates are unaf- 
fected. For an IH, the situation is reversed. Once again, these results are apparent in 
the plots shown in Figure |21 

In a charge discriminating detector, which detects /i~ and fi'^ rates individually, a com- 
parison of the fi~ event rate with that expected from vacuum oscillations is enough to de- 
termine sign(A3i), whether the signal is measured in the ranges 1 or in ranges 2 28 . Water 
Cerenkov detectors are insensitive to the sign of the leptonic charge on an event-by-event ba- 
sis. Therefore, we sum over the rates of /i~ and /i"*" events and label the sum as 'muon events'. 
Hence, for both normal and inverted hierarchies, the total muon event rate due to matter 
oscillations will be less(more) than the corresponding rate in the case of vacuum oscillations 
in the case of ranges 1(2). 

While the summation of events over muons and anti-muons dilutes the sensitivity of 
water Cerenkov detectors to matter effects compared to charge discriminating detectors, the 
proposed megaton mass overrides this disadvantage and provides the statistics necessary for a 
determination of the hierarchy. It is also to be noted that the /i~ event rates are 2 to 3 times 
larger than those of the /i"*" rates due to the larger u^ — N cross section. Thus a larger change 
in the total muon event rate and a statistically stronger signal of matter resonance effects is 
envisioned in the case of a NH than in the case of an IH in both types of detectors. 

The total number of muon (or anti-muon) charged current (CC) events can be obtained 
by folding the relevant CC cross section with the survival probability, the incident neutrino 
fiux ,the efficiency for detection, the detector mass and the exposure time. In our calculations. 
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the cross sections for i/^ — N and z/^ — N interactions have been taken from llHj. The total CC 
cross section is sum of quasi-elastic, single meson production and deep inelastic cross sections. 

Before describing our calculation and the results, we first state the assumptions we have 
made. We assume that the neutrino oscillation parameters IA31I and ^23 will be determined 
to better than 10% accuracy by the long baseline experiments MINOS ^Hj and T2K [S] by 
the time megaton size water Cerenkov detectors are operative. In addition, we have assumed 
that similar precision will exist for the measured values of A21 and 612 EHI , even though 
the dependence of our results on them is marginal. Based as they are on survival rates and 
large matter effects, our conclusions are also effectively independent of the CP violating phase 
6c p- Finally, for our method to yield results over running times of 3-4 years, sin^ 26'i3 must 
be > 0.05. 

For the first set of our calculations, we use the following values for all neutrino param- 
eters known to be non-zero: IA31I = 0.002 eV2, A21 = 8.3 x lO'^ eV^, sin^ ^23 = 0.5 and 
sin^6'i2 = 0.28. We calculate the event rates for four different values of 6*13, i.e. sin^ 26*13 = 
0.05,0.1,0.15,0.2. In a later set of calculations, we vary 6*23 and IA31I within the ranges al- 
lowed by SK atmospheric neutrino data and K2K results to gauge the dependence of our 
results on them. 

To study the sensitivity of the megaton water Cerenkov detectors to the hierarchy, we 
compute the expected number of /i~ + /z"*" events in the following four cases: 



• No oscillations, Ndowii, 

• Vacuum oscillations, Nvac, 

• Matter oscillations with NH, N^h and 

• Matter oscillations with IH, Nm. 



If the difference between Nnh and Njh is statistically significant, then sign(A3i) can be assumed 
to established with this same significance. Thus 

|Nnh-Nih| 
= /N— 

gives the sigma confidence level with which the expectation for NH differs from that for IH. 

We note that Eq.© tacitly assumes that atmospheric neutrino fluxes will be measured 
well enough in the next decade to allow use of the absolute event rates. However, if this is 
not the case, one can also use the measured number of downward going events, which have 
the same corresponding values of | cos6'zenith| as the upgoing events, as an estimate for the 
number of events in case of no oscillations. This is the reason why we have denoted these 
as Ndowh above. In such a situation, we can usefully consider the ratios Rnh = NNH/Noown 
and RiH = Nih/Ndowii- The statistical significance of the difference between the ratios is then 
defined as 

N., = \^-^-^^^\ . (6) 

R-NHV 1/NnH + l/Noown 
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Table 1: Integrated muon event numbers for SK (337.5 Ktyr) for the energy range E = 5 - 10 GeV and 
pathlength range L — 6000 - 9000 km. The expected event numbers in the case of vacuum oscillations (Nvac ) 
and matter oscillations for both NH (Nnh) and IH (Nm) are calculated with sin^ 023 — 0.5, IA31I = 0.002 eV^, 
and various different values of 6*13. The number of downward events Ndowh = 178 is used in calculating the 
ratios Rnh = Nnh/Ndowii and Rm = Nm/NDown 



Our first set of results is for SuperKamiokande (SK) for 15 years of running. Since the 
SK fiducial volume is 22.5 Kt, this gives us an exposure^ of 337.5 Ktyr. We use the Bartol 
fluxes from jJU] and incorporate the SK cross sections, resolutions [lUll^ and efficiencies |52j . 
In calculating the number of events, we have integrated over the neutrino energy range = 5-10 
GeV and pathlength range L = 6000 - 9000 km. For A31 = 0.002 eV^, the resonance in earth's 
mantle occurs for energy of about 6.5 GeV. Since we are looking for resonant amplification 
of matter effects, we have chosen the energy range so that it encompasses this energy. We 
note that this choice is also dependant on the fact that IA31I = 0.002 eV^. For other values of 
IA31I, the appropriate energy range will be different. It was shown earlier that the pathlength 
for which the resonant matter effects are maximum is about 7000 km. In this calculation, we 
choose a range of pathlengths, 6000 - 9000 km, in order to utilize the enhancement associated 
with this value. 

Table^gives our results for SK for an exposure time of 337.5 Ktyr. We find that, even for 
the largest allowed value of sin^ 2^13, the expectations for NH and IH differ by 2.8cr if we have 
reliable predictions for neutrino fluxes. If one needs to take ratios to cancel flux uncertainties, 
then the confidence level in the difference reduces to 2.2a. These low confidence levels result 
from the small event rates of atmospheric neutrinos at high {i.e several GeV) energies. Since 
the observation of the resonant amplification of matter effects requires energies in the 4 — 10 
GeV range, the larger exposures which are possible at planned megaton size detectors such as 
HyperKamiokande (HK) ESI (higher by a factor of 6 to 20 compared to SK) are required 
to obtain a statistically significant signal. The proposed fiducial volume of HK is 545 Kt. 
Hence, with a running time of only 3.3 years, it will be possible to have a total exposure of 1.8 
Mtyr, which is a factor 6 larger than what is possible at SK. In the calculation of the muon 
event rates for HK, we have used the efficiencies and resolutions of SK 

Figure El shows the L/E distribution of events for matter oscillations (for both a NH and 
an IH) by plotting the events versus Logio(L/E) for HK in the energy range from 5-10 
GeV and pathlength range 6000 - 9000 km, with sin^ 2^13 = 0.1. Assuming an L/E resolution 
similar to SK allows a division into five bins in Logio(L/E) for this L and E range [SI]. The 
drop in Nnh relative to Nm is evident in the third and fourth bins in this figure. In these 

^In the latest SK paper the exposure is 92 Ktyr corresponding to 1489 days live time from May 1996 
to July 2001. Preliminary results for SK-II include a further 627 days live time [13 ■ 
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Table 2: Integrated muon event numbers for HK (exposure of 1.8 Mtyr) for the energy range E = 5 - 
10 GeV and pathlength range L = 6000 - 9000 km. The expected event numbers in the case of vacuum 
oscillations (Nvac) and matter oscillations for both NH (Nnh) and IH (Nm) are calculated with sin^ 023 = 0.5, 
IA31I = 0.002 eV^, and various different values of 6113. The number of downward events Noown = 951 is used 
in calculating the ratios Rnh = Nnh/Ndowh and Rm = Nih/Ndowh 
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Figure 3: Muon {i 



event distribution for HK(1.8 Mtyr) of matter oscillations for NH and for IH 



plotted vs Logio(L/E) for the E and L range E = 5-10 GeV, 
used are sin^ 6*12 = 0.28, sin^ 6*23 = 0.5, sin^ 26*13 = 0.1, IA31I = 



L = 6000-9000 km. The values of parameters 
2 X 10-3 eV^ and lAail = 8.3 x 10"^ eV^. 



energy and zenith angle ranges the change in P^^ is predominantly due to the matter effect in 
P^e- In Table El we present the N^- values calculated using EqlHl We first sum over the events 
of all the 5 bins shown in Figure El and then compute No-. For sin^ 2^13 = 0.1 the sensitivity 
is 3.7(7. If instead we consider only the 3rd and 4th bins of Figure |21 and sum the number 
of events in these two bins then we get Nnh = 400 and Njh = 475, which correspond to a 
sensitivity of 3.8(7. Making a narrow choice of bins in (L/E) does not improve the sensitivity 
as the effect is distributed over a larger L and E range. However, we will show later that for 
the L and E ranges for which the matter effects in P^^ arise dominantly from P^^, binning in 
(L/E) leads to a great improvement in the sensitivity to mass hierarchy. 

As shown in Tabled we obtain a ^ 3.7(7 difference between Nnh and Nih for sin^ 2^13 ^ 0.1, 
for an exposure of 1.8 Mtyr. Even when we have to take ratios with respect to downward 
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0.5 
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490 571 3.7a 


0.52 0.6 2.8a 


0.6 


596 


475 576 4.6 a 


0.5 0.6 3.6 a 



Table 3: Integrated muon event numbers for HK (exposure of 1.8 Mtyr) for the energy range E = 5 - 
10 GeV and pathlength range L = 6000 - 9000 km. The expected event numbers in the case of vacuum 
oscillations (Nvac) and matter oscillations for both NH (Nnh) and IH (Nm) are calculated with sin^ 20i3 = 0.1, 
IA31I — 0.002 eV^, and various different values of 6*23. The number of downward events Noown = 951 is used 
in calculating the ratios Rnh = Nnh/Ndowh and Rm = Nih/Ndowh 



events, this difference is close to 3a. With an exposure of about 3 Mtyr (which corresponds 
to a running time of about 6 years), one can obtain a 3a difference between Nnh and Nih, for 
sin^ 2^13 ~ 0.05. 

Next, we check how N^- and N^-p^ change with variation of neutrino parameters. We first 
illustrate the dependence on 6*23. For this purpose we fix sin^2^^i3 = 0.1 and consider three 
different values of sin^ ^23 = 0.4, 0.5 and 0.6. The sensitivities for these values of ^23 are shown 
in Table El for HK exposure of 1.8 Mtyr. The energy range (5 - 10 GeV) and the pathlength 
range (6000 - 9000 km) remain the same because IA31I is unchanged. We note from Table El 
that the sensitivities are better for sin^ 6'23 = 0.6 and worse for sin^ 6'23 = 0.4 compared to the 
sin^ 6'23 = 0.5 case. The reason for this is simple. As mentioned earlier, the matter effects in 
P^^ arise mostly due to the matter effects in P^e- The matter term in P^e is proportional to 
sin^^23- Therefore, for a larger value of sin^ ^23 we get a larger change in P^^ due to matter 
effects and hence there is a larger difference between the expected rates for NH and IH cases. 

We next study the change in the sensitivities with variation in IA31I. We fix sin^ 6'23 = 0.5, 
sin^ 26113 = 0.1 and consider three different values of IA31I = 0.001,0.002 and 0.003 eV^and 
calculate the muon event rates for HK exposure of 1.8 Mtyr by integrating over appropriate 
energy and pathlength ranges. From Eq.(0), we note that the resonance occurs at different 
energies for different values of IA31I. For 0.003 eV^, it occurs at 9.5 GeV. So we choose the 
energy range of integration to be 7— 12 GeV. For 0.001 eV^, the resonance occurs at 3.2 GeV, 
so we choose the energy range to be 2 — 5 GeV. But the pathlength range remains unchanged 
at 6000 - 9000 km. This is because Eq.(j2I) shows that the expression for the pathlength for 
resonant amplification of matter effects is independent of A31. The event numbers and the 
corresponding sensitivities for this case are given in Table IH We note that in general one 
expects a fall the in sensitivity for higher values of IA31I due to the higher resonance energies 
(and consequently lower number of total events) ^ . 

It was shown in references |40| l^ that for certain ranges of L and E, P^^ is larger compared 
to its vacuum value, due to a sharp fall (by as much as 70%) in P^^. For IA31I = 0.002 
eV^, this occurs in the energy range 4-6 GeV and pathlength range 8000 - 10500 km. In 
Figure m the total muon event rates are computed and plotted as function of L/E for values 

^That this is not the case for IA31I = 0.002 compared to IA31I = 0.001 in our table is attributable to the 
fact that the energy range for the best fit value of IA31I — 0.002 has been optimized, while that for the other 
two values has been chosen merely for illustrative purposes. 
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Table 4: Muon event numbers for HK for an exposure of 1.8 Mtyr for different values of | A31 1. The pathlength 
range, L — 6000 - 9000 km, is the same for all three cases but the energy ranges are different because the 
resonance energies are different. The expected event numbers in the case of vacuum oscillations (Nvac) and 
matter oscillations for both NH (Nnh) and IH (Nih) are calculated with sin^ 6123 — 0.5 and sin'^ 2^13 — 0.1 
Because of different energy ranges, the number of downward events in each case is different and are listed in 
table. They are used in calculating the ratios Rnh = Nnh/Ndowh and Rm = Nih/Ndowh 



sin^ 2^13 


N 

vac 


Nnh Nih N, 


Rnh R-ih Ncrp^ 


0.05 


40 


68 46 2.6 a 


0.21 0.14 2.5 a 


0.10 


41 


96 58 3.9 a 


0.3 0.18 3.4(T 


0.15 


49 


130 74 4.9 a 


0.4 0.23 4.1 a 


0.20 


50 


158 86 5.7 a 


0.49 0.26 4.8 a 



Table 5: Integrated muon event numbers for HK (exposure of 1.8 Mtyr). The events selected are those in 
the energy range E = 4 - 8 GeV and pathlength range L = 8000 - 10500 km, with Logio(L/E) = 3.2 - 3.3. 
The expected event numbers in the case of vacuum oscillations (Nvac) and matter oscillations for both NH 
(Nnh) and IH (Nm) are calculated with sin^ 6*23 — 0.5, IA31I = 0.002 eV^, and various different values of 
013 . The number of downward events Noown = 325 is used in calculating the ratios Rnh = Nnh/Ndowh and 

RlH = Nih/Ndowii 



of Logio(L/E) in the range 3.0 - 3.4. If we assume SK resolutions, this range can be divided 
into four bins. We see that for the bin with Logio(L/E) = 3.2 — 3.3 the change due to 
the sign of A31 is particularly pronounced. Note also that, as mentioned above, unlike the 
previous ranges, here the expected number of muon events for the case of NH is larger than 
that for the case of IH. This phenomenon, however, occurs only for narrow ranges in energy 
and pathlength and hence one needs a detector with good L/E resolution to observe it. Given 
a detector with such resolution however, the difference in the expected muon rates in the bin 
with Logio(L/E) = 3.2 — 3.3 is an excellent indicator of the neutrino mass hierarchy. From 
the event numbers given in Table El we see that the difference in this bin corresponds to a 4a 
signal for sin^ 2^13 = 0.1 with an exposure of just 1.8 Mtyr in HK. If we sum the events in the 
two bins with Logio(L/E) = 3.1 — 3.2 and 3.2 — 3.3, the difference between the expectations 
for NH and IH corresponds to a sensitivity of 3a. 

Table El also gives the number of events in matter (A31 > 0) and the values of a sensitivity 
for this Logio(L/E) bin for three values of 613. This shows a higher sensitivity as compared to 
Tableland more than 2.5cr sensitivity can be achieved for sin^ 2^13 = 0.05 with the exposure 
considered here. 
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Figure 4: Muon (i^^ + P^) event distribution for HK(1.8 Mtyr) in matter and in vacuum plotted vs Logio(L/E) 
for the E and L range E = 4-8 GeV, L = 8000-10500 km. The values of parameters used are sin^ 612 = 0.28, 
sin^eza = 0.5, sin^ 26'i3 = 0.1, IA31I = 2 x lO^^ eV^ and IA21I = 8.3 x lO^^ eV^. 



4 Conclusions 



• The neutrino mass hierarchy, presently unknown, is a powerful discriminator among 
various classes of unification theories. We have shown that the total muon (z/^ -|- z/^) 
survival rate in atmospheric events can provide a novel and useful method of determining 
the hierarchy in megaton water Cerenkov detectors^. 

• The rates for both normal and inverted hierarchies differ from the vacuum rate and 
from each other due to large matter sensitive variations not only in P^c but also in P^^. 
We have identified energy and pathlength ranges where these effects can be fruitfully 
observed in a statistically significant manner for running times of ~ 3 — 4 years for values 
of sin^ 013 > 0.05. 

• We have calculated the a sensitivities for normal and inverted hierarchies via differences 
in absolute events as well as differences in event ratios (to eliminate uncertainties in the 
atmospheric flux) for various choices of the neutrino parameters. We find that the energy 
ranges 5 — 10 Gev and 4 — 8 GeV, when combined with the length ranges 6000 — 9000 km 
and 8000 — 10500 km respectively^^, allow for a statistically significant determination 
of the hierarchy. The first (L,E) range combination allows us to probe the P^e-induced 

^Even though we have presented the results for HK, our conclusions should remain valid for similar detectors 
such as UNO or the one proposed in Frejus tunnel 34 in Europe. 

^°In the latter case, as discussed in the context of Figure 01 the effect is apparent in the L/E bin specified 
by Logio(L/E) = 3.2 - 3.3. 
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matter effects in P^^, while the second combination provides a way of detecting the P^^ 
effects in P^^. 

• Finally, we have also studied the dependence of our results on the uncertainties in the 
presently known values of sin^ 6'23 and |A3i|. 
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